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Abstract 

We discuss here the possibility of explaining the rotational velocity curves of the 
galaxies in f{TZ) gravity by solving field equations numerically. In a given constant 
rotational velocity region, we prove that all values of rotational velocities does not lead 
to an analytic solution of the field equations in vacuum. We then obtain the numerical 
solutions of the field equations which suggests that a very slightly modification from 
linear relations of TZ can be an alternative to the so called dark matter. 
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1 Introduction 

There is no priori reason to consider that the gravity is linear in Ricci scalar TZ except 
it best matches the observation (excluding dark matter and dark energy). Dark energy is 
related to the mystery of accelerated expanding Universe whereas dark matter is related 
to the observed rotation velocity curves of the galaxy. Here, we address the problem of 
the dark matter. The present Newtonian dynamics and Einstein’s gravity can not explain 
these phenomenas which leads to idea that either one or both sides of the Einstein held 
equations (EFE) are incomplete. The right hand side of the EFE comprise of the matter 
part whereas left side is the geometric part. Einstein held equations simply tells that if we 
put some matter in the space, the space counter acts to tell the matter how to react in a 
given geometry. There are many proposals that address the dark matter problem. A few 
of them are weakly interacting massive particles (WIMP) [HI E], scalar held from extensions 
of SM [7], modihed Newtonian dynamics (MOND) [10], modihed gravity [161 ESI ll]- The 
modihed gravity modihes the left hand side (geometric part) of the Einstein held equations 
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(field equations then called). It then can be used to explain the observed rotational velocity 
curves to be a geometric effect only. 

The mysterious dark matter is the type of matter supposed to be surrounding galaxies 
that give the observed rotational velocity curves. Einstein Hilbert action could also be one of 
the reason that there is some discrepancy between observed and expected tangential velocity 
of a particle moving around the galaxy. 

In this article, we present a simple version of /(7^) gravity model by solving held equa¬ 
tions numerically that can describe the dark matter as a geometric effect only. The scheme of 
the article is as follows; using a metric that can describe the constant velocity region of space 
around a galaxy we try to solve the held equations for the metric in f{TZ) gravity analyti¬ 
cally and numerical results of the held equations suggest that very slight modihcation from 
Einstein’s gravity could explain the discrepancy between observed and expected rotational 
velocity curves. The article is as follows; hrst a metric that can describe the constant velocity 
region of space around a galaxy is given. We then move on to solve the held equations for 
the metric in f{TZ) gravity analytically. The numerical results of the held equations succeed 
that. 


2 fi'TV) gravity field equations in vacuum 

The metric which describes the static and spherically symmetry in the constant rotational 
velocity regions is 

+ r^sin (1) 

The metric coefficient z/(r) can be calculated from the Euler Lagrange equations of the 
Lagrangian of the metric. From Euler Lagrange equations, we get i>{r) = 2mln(r/ro) where 
m = Vtg being the rotational velocity of a particle around galaxy, c is the speed of 

light and tq is the constant of integration obtained in deriving 
The /(7^) modihed gravity action is written as 

s = / /(K) dV , (2) 

here /(7^) is an arbitrary analytic function of Ricci scalar TZ. The variation of the above 
action with respect to the metric gives the following held equations 

F(7^)7^^, - - (V^V, - g^,D) F{R) = 0 , (3) 

where F(Jl) = df /dTZ. The contraction of the above equation gives 

F(7^)7^ - 2/(7^) + 3nF{R) = 0 . (4) 

Upon using eq. (|^ in eq. ([^ we get modihed held equations as HU 

F{lZ)R.gy - VgS/yF = ^gg^ (F7^ - DF ), (5) 
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( 6 ) 


also differentiation of eq. Q with respect to ‘r’ gives 

KF' - 7^'F + 3 (DF)' = 0 . 

Where ' represents derivative with respect to ‘r’. If we use modihed held equations to hnd the 
solution then the solution must also satisfy eq. Q to be a solution of original held equations 
given by eq. (|^. From eq. ([^, we see that — V^) /gfj,^ is independent of 

the index /r, thus = 0. From the relation A^ — Ay = Q we can write three equations 

(Aq — Ai = 0, Ai — A 2 = 0 and Aq — A 2 = 0). The equations Aq — Ai = 0, Ai — A 2 = 0, 
the rr component of the held equations and eq. Q can be written as 

F" - hu + A')F' + -(z/' + A')F = 0 , (7) 

2 r 


/2 1 ^ , 2A / , 2 „F" 2\ F' 


Z2 +u -;;:(z2+-)(z/+A)- — + 2— - 2— - ( A + - ) — , 


r ) F 


f = Fe 


-A 


U 


1 w ' f ' 4A F" 

2" A+^>-27 + (‘'+; W 


( 8 ) 

( 9 ) 


f 

n = 2 ^ -3e"^ 
F 


F fl / , 2A F 

—+ (-b -^) + ;) -7 


( 10 ) 


respectively. 

Introducing rj = ln(r/r*) with r* = 1, eq. 0-0 now become 


d^F 


1 /dz/ dAAA dF /dz/ dA' 

2 \dr] df])) df] \dri dr]/ 


( 11 ) 
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( 12 ) 


/ = Fe-^-2^ 


dV _ dz/ _ 1 /dz/ dAA dz/ _ ^dA /dz/ A 1 dF 

dr;^ dr; 2 \ dr; dr; / dr; dr; \ dr; / F dr; 


(13) 


TZ = 2j- 


1 d^F _ 1 dF /I /dz/ _ dAA A 1 dF' 
F dr;^ F dr; \2 \dr; dr;/ "*” / F dr; 


(14) 
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With the introduction IjF {dF/dvi) = m, eqs. (11) and (12) simplihes to 


dti 


1 /dz/ dA\ 




du dA 


( 15 ) 


d^_(^ /^y_i/^ 2^— —^ 

dp^ dr; l^dr ;) 2 l^dr; ) l^dr; dr; J 


dA , , 


( 16 ) 


even though above equations are second order differential equations in u but since u has 


already been given as a function of r earlier, the introduction of u makes eqs. (11) and (12) 
hrst order with the increase in number of equations. Using du/drj + from eq. (16) into 


eq. (15) we get 


d^v _dv /^dz/y_/ ldv\fdv dA 
dr;^ dr; \ dr; y \ 2 dr; J \ dr; dr; 


/ \\ / 1 dz/A /dz/ dAA 

2(l-e)+2(l---j« + 2(- + -j=0, 


(17) 


Eq. (17) and eq. (15) or eq. (16) provide the solution of A and u (z/ is known). 


3 Dark matter as a geometric effect 


To obtain the geometric interpretation of the constant velocity regions, we solve eq. (15|) 
and eq. 0 for A and u with u = 2m{r] — rjo). Eq. ( [I5| ) and eq. 0 can now be written 


as 


drr , 


^ IdAA , dA ^ 

1 + m + -— ] u + 2m + — = 0 , 
2 dr] J dr] 


( 18 ) 


—2m + dm'^ — (1 + m) ^2m + — 2 (l — e^) + 2 (1 — m) rr + 2 ^2m + = 0 . (19) 


Using u and du/dr] from eq. (19) into eq. (18) we obtain 

Id^A 1 /3 , AdA 1 /dAA 

- e — m(l — 'fn) -5 1 


2 dr;^ 1 — m ^2 " J dr] 2 \dr] J 

(1 — e^{l + m^) — 1 


( 20 ) 


+ 2m = 0 , 


(1 — m)^ (1 — m)^ 

This is a second order differential equation in A. After obtaining A from this equation 


one can use either eq. (18) or eq. (19) to obtain u. The obtained A and u must satisfy 
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the third equation (which has not been used in getting A and u). As discussed in [T], the 
tangential velocity of a test particle in stable circular orbit around the galactic center is 
about 200 — 300Km/s [HI [121 E] thus m k, 0(10“®). Hence, we must get an approximate 
function of Ricci scalar around m ~ 0(10“®). 

Introducing new variable as dA/dr/ = p then we can write d/dp = p d/dA. This allows 
us to write eq. (20) as 


2 dA 1 


m 


2 ® 


m(l — m) \p - p 


(1 — e^Y e^{l + rrr) — l 


( 21 ) 


(1 — mY 


(1 — mY 


+ 2m = 0 


Now, using the transformation p = l/uj and = P, eq. (21) after simplihcation can be 
written as 


do; 

1 

to 

1 

isD 

1 
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This is the nonlinear hrst order Abel differential of hrst kind of the form 

^ = p{9)u^ + q{9)u'^ + r{9)u + s{9) , 

du 


with 


m 2(1 — 9Y 2(6*(l + m^) —1) 
P(^) = ^ 


q{9) = 3 - 
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2m(l — m) 
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If we choose that the particular solution pp of the above equation satishes 
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then eq. (22) becomes Abel differential equation of second kind [S], of the form 
m 2(1 — 9Y 2(6*(1 + m?) — 1) 


da; 


_9 9(1-mY 
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In the case when. 


pW 




= - 3p(«)r(») , 


( 22 ) 


(23) 


(24) 


(25) 
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where is a arbitrary constant then the general solution of eq. (22) is given by [H], 


p[e) 


-q{e) ± ^Jq{9Y - ?>p{e)r{d) 


3p(0) 


( 26 ) 


where V satishes the separable differential equation 


dl/ 

~de 


(v^ + v^ + Kv) 




pW 


(27) 


To get the solution, we must satisfy the condition given by eq. (25). The plot of 

P{9) 


1 


a = 


^q^ey -‘ip{e)r{e) d^ \ ^jq{ey -^p{e)r{e) 


is given in £g. ([^ and fig. (|^. The plot allows us to understand when (i.e. for what m and 
6) could we get an analytic solution to the eq. (22). 
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Figure 1: When 6^ = 0 we have eq. 
arbitrary negative large value of A(r). 


(25) satished 


for all values of m. 


This demands an 
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Figure 2: The plot is plotted for all possible values of 9 (i.e. 0 < 6^ < cxd). When we have 
m = 1 eq. (25) is obeyed for all values of 6 except —oo. For all other values of m if we take 


the domain of 6* to be 0 < 6* < 1 (i.e. all possible values, (0, oo) is being mapped onto (0,1) 
with the mapping function then we observe the violation of eq. (25). 


It has also be proven in [T5] that A = constant is not the solution of eq. (18) and (19) in 
the limit that and higher order terms in m can be neglected. 


4 Numerical solution to the field equations 


From the hgures ([^ and (|^, any analytic solution of eq. (22) in the presented scheme around 
m = 0(10“®) can not be obtained for all possible values of 6. We then hnd the numerical 
solution to the held eqs. ([^ and ([^. The initial conditions taken were 


= 10 Kpc , F{rini) = 1 , 


F'i 


= 0 


A(rini) = 2 X 10 


-6 


The initial conditions were chosen such that initially our model of fiTZ) gravity matches 
onto Einstein’s gravity (the value of A was chosen by brute force attack). 
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Figure 3: The radial coefficient of the metric. As we move away from the center of the galaxy 
A(r) decreases till r 300 Kpc then it starts to increase. 

From the fig. (|^ we see that the metric coefficient A(r) decreases initially till r 
300 Kpc. For r > 300 Kpc we see a monotonically increasing function. 



r (Kpc) 

Figure 4: Normalized Ricci scalar as a function of the radial parameter r, moving away from 
the center of the galaxy we have 7^ —)■ 0. 

The Ricci scalar TZ decreases as we move away from galactic center. As we get around 
r 100 — 200 Kpc, Ricci scalar goes to zero showing that we have reached the galactic edge. 
This should be the case as away from galaxy the volume of a geodesic sphere should match 
on to that of a ball in Euclidean space. 
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r (Kpc) 

Figure 5: Normalized / as a function of the radial parameter r, moving away from the center 
of the galaxy we have / —)• 0. 

The / as a function of r is given in fig. ([^. In moving away from the galactic center / 
decreases and goes to zero thus matches with IZ = Q case. 



Figure 6: The blue thick dashed line is our model whereas the solid black line is f{TZ) = TZ. 
Moving slightly away from our f{TZ) model starts to deviate from linear relations very 
slightly (a clearer version of this statement is given in table Q). As we reach the boundary 
of the galaxy we have 7?. ~ 0 and f{TZ) ~ 0. 

The fig. (|^ is the parametric graph of f{TZ) and TZ. We see from the graph as we 
move away from the galaxy the line of f{TZ) is approximately linear. The solid black line 
corresponds to f{TZ) = TZ, it overlaps our model. The black line has been drawn to compare 
Einstein’s gravity and our f{TZ) gravity model. Due to very slight modihcation from linear 
relation it is very hard to see any difference between f{TZ) = TZ and our model. But we give 
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some numerical values to see that difference. 


r (Kpc) 

Normalized Abs [77] 

Normalized Abs[/] 

10 

1 

1 

12 

0.694395 

0.694347 

15 

0.444373 

0.444302 

20 

0.249931 

0.249868 

40 

0.0624127 

0.0623267 

70 

0.0203219 

0.0202357 

100 

0.00991653 

0.00983378 

200 

0.00241698 

0.00233415 

500 

0.000314893 

0.000229789 


Table 1: Some numerical values of TZ and /(7^) in our model. The model starts from 
f{TZ) = TZ, moving away from galactic center it deviates from Einstein’s gravity. ‘Abs’ 
represents the absolute value. The precision of the numerical values is the machine precision 
of the MATHEMATICA software set at 15.9546 whereas the accuracy goal is set to 10. Thus, 
we expect the numerical error in each result to be less than or equal to 10“^° + 
where |x| here is the TZ at different values of r. 

From the table Q, a very slight deviation from Einstein’s gravity can explain the constant 
rotational velocity of a particle in moving away from the galaxy. The table [T] suggests that 
the a very slight negative correction to the f(JZ) = 7Z can explain the constant galactic 
rotation curves i.e. fiJZ) < IZ. 

5 Conclusions and discussion 

In this article, we started of with a metric that describes the constant rotational velocity 
region in spherically symmetric and static space-time. We then obtained the numerical 
solution to the held equation to get an appropriate estimate of f(JZ) gravity that could be an 
alternative to the dark matter. The resultant graphs showed that a very slight modification 
to the Einstein’s gravity can be an alternative to the dark matter. As said in the last 
section, the correction up to r = 500Kpc we find is negative i.e. f(JZ) < 7Z, shown in 
table Although we are not giving any numerical result beyond r = 500Kpc but it is 
worth mentioning that in our numerical solutions after r ^ 885Kpc, the modifications in the 
f{TZ) = TZ becomes positive i.e. for r ^ 885Kpc, f{TZ) > TZ. 

The weak field limit of the f{TZ) gravity models has been discussed for star like objects 
in [31 Eli. If we assume that f{TZ) is an analytical function at the constant curvature Rq, 
that rrifjj/r 1, where is the effective mass of the scalar degree of freedom of the theory, 
and that the fluid is pressureless, the post-Newtonian potentials ^(r) and <I)(r) are obtained 
for a metric 

dS'^ = — (1 — 24/(r)) + (1 + 2<I>(r)) dr^ -|- r^dfl^. 
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In that way, the behavior of \l/(r) and <I)(r) outside the star can be estimated. This then 
gives the value of the post-Newtonian parameter 7 to be 1/2 [1]. From the solar system 
observations, we know that 7 = 1 [I]. This inconsistency between expected and measured 
value of 7 then is used to rule out most of the modihed gravity models. 

T. Chiba, T.L. Smith and A.L. Erickcek have done the same analysis for the /(7^) gravity 
function of the form f(Jl) = {7l/aY~^^ in [3]. They concluded that in the their section III. 
Case Studies “...this analysis is incapable of determining whether fiJZ) = gravity with 
6^1 conflicts with solar system tests.”. Thus, we can conclude that results presented in 
fig- (|^ and table Q does not contradict with the solar model tests. 

One should keep in mind that the graphs presented above are obtained in approximation 
that the rotational velocity of a particle around the galaxy in radial direction remains con¬ 
stant however this is not exactly the case observed but much close to the observations, trend 
of rotational velocity curve for different galaxies is different but slightly. The assumption is 
valid for at least a signihcant region of the total tangential velocity prohle. If one can estab¬ 
lish an exact relation between the rotational velocity and distant from the galactic center 
then one can obtain more accurate numerical results of the held equations. In that case, we 
might even also be able to hnd the exact solution to the held equations. 

In conclusion, we found that regions with hat galactic rotations curves do not require any 
kind of dark matter if we try to explain these in f{TZ) gravity. The constant rotation curves 
are then a consequence of the additional geometrical structure provided by the modihed 
gravity function /(T?.). We found that a very slight modihcation to the Einstein-Hilbert 
Lagrangian may account for the existence of “dark matter”. In the fiTZ) gravity modihed 
theory in principle we can rewrite the held equations in terms of the Einstein tensor TZ and 
interpret the remaining term(s) as a geometrical energy-momentum tensor which then gives 
us the option not to include the exotic type of matter called dark matter. In this way, the 
presence of the remaining terms can provide us with an elegant geometric interpretation of 
the dark matter. 


6 Acknowledgements 

This work is supported by National University of Sciences and Technology (NUST), Sector 
H-12 Islamabad 44000, Pakistan. 


References 

[ 1 ] Christian G. Bohmer, Tiberiu Harko, and Francisco S.N. Lobo. Dark matter as a 
geometric ehect in gravity. Astroparticle Physics, 29(6):386 - 392, 2008. 

[ 2 ] Annamaria Borriello and Paolo Salucci. The dark matter distribution in disc galaxies. 
Monthly Notices of the Royal Astronomical Society, 323(2):285-292, 2001. 


11 


[3] Takeshi Chiba, Tristan L. Smith, and Adrienne L. Erickcek. Solar system constraints 
to general /(r) gravity. Phys. Rev. D, 75:124014, Jun 2007. 

[4] Takeshi Chiba, Tristan L. Smith, and Adrienne L. Erickcek. Reply to “comment on 
‘solar system constraints to general /(r) gravity’”. Phys. Rev. D, 77:108502, May 2008. 

[5] Valerio Faraoni and Nicolas Lanahan-Tremblay. Comment on “solar system constraints 
to general /(r) gravity”. Phys. Rev. D, 77:108501, May 2008. 

[6] Rudy C. Gilmore. Mass limits on neutralino dark matter. Phys. Rev. D, 76:043520, 
Aug 2007. 

[7] I. F. Ginzburg, K. A. Kanishev, M. Krawczyk, and D. Sokolowska. Evolution of the 
universe to the present inert phase. Phys. Rev. D, 82:123533, Dec 2010. 

[8] Gerard Jungman, Marc Kamionkowski, and Kim Griest. Supersymmetric dark matter. 
Physies Reports, 267(56): 195 - 373, 1996. 

[9] M.K. Mak and T. Harko. New method for generating general solution of abel differential 
equation. Computers & Mathematies with Applications, 43(12):91 - 94, 2002. 

[10] M. Milgrom. A Modihcation of the Newtonian dynamics as a possible alternative to the 
hidden mass hypothesis. Astrophys.J., 270:365-370, 1983. 

[11] T. Multamaki and I. Vilja. Spherically symmetric solutions of modihed held equations 
in /(77) theories of gravity. Phys. Rev. D, 74:064022, Sep 2006. 

[12] Massimo Persic, Paolo Salucci, and Fulvio Stel. The universal rotation curve of spiral 
galaxies i. the dark matter connection. Monthly Notices of the Royal Astronomical 
Society, 281(l):27-47, 1996. 

[13] R. Saffari and Y. Sobouti. An f(TZ) gravitation for galactic environments. A&A, 
472(3):833, 2007. 

[14] P. Salucci, A. Lapi, C. Tonini, G. Gentile, I. Yegorova, and U. Klein. The universal 
rotation curve of spiral galaxies ii. the dark matter distribution out to the virial radius. 
Monthly Notices of the Royal Astronomical Society, 378(1):41-47, 2007. 

[15] Muhammad Usman. Note on “dark matter as a geometric effect in f{TZ) gravity”. 

[16] Y. Sobouti. An /(77) gravitation for galactic environments. A&A, 464(3):921-925, 2007. 


12 



